Engineering multiferroism in CaMnOs 
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Prom first-principles calculations, we investigate tlie structural instabilities of CaMnOs. We 
point out that, on top of a strong antiferrodistortive instability responsible for its orthorhombic 
ground-state, the cubic perovskite structure of CaMnOs also exhibit a weak ferroelectric instability. 
Although ferroelectricity is suppressed by antiferrodistortive oxygen motions, we show that it can 
be favored using strain or chemical engineering in order to make CaMnOs multiferroic. We finally 
highlight that the FE instability of CaMnOs is Mn-dominated. This illustrates that, contrary to 
the common believe, ferroelectricity and magnetism are not necessarily exclusive but can be driven 
by the same cation. 



The family of ABO3 perovskite oxide compounds con- 
situtes an important class of multifunctional materials 
which, within the same simple cubic structure, exhibit 
a wide variety of behaviors and range from ferroelectric, 
piezoelectric or non-polar anti-ferrodistorted insulators 
to metals or superconductors, eventually combining an 
additional magnetic order p]. Recently, the search for 
magneto-electric multiferroics operating at room temper- 
ature has motivated the study of oxides combining ferro- 
electric and magnetic orders like BiFeOa or YMnOs j2ll3]. 
In this framework, the case of CaMnOa, a magnetic in- 
sulator which remains paraelectric, was considered as a 
prototypical example to reveal the antagonist role of par- 
tial d-shell occupacy in respectively favoring magnetism 
but penalizing ferroelectricity j^. 

CaMnOa is a G-type antiferromagnetic insulator[5]- It 
cristallizes at low temperature in a paraelectric Pnma or- 
thorhombic ground-state structure which can be viewed 
as a small distortion of the cubic perovskite structure, 
produced by the tilting of oxygen octahedra (a~b+a~)|6] 
[3. Because of the Mn''"'" configuration, CaMnOa does not 
exhibit Jahn- Teller distortion : in the cubic structure, 
the 5d orbitals of Mn ion split into 3-fold degenerated 
t2g levels and two-fold degenerated e^ levels with occupa- 
tion being tjg and e^p . CaMnOs does not either exhibit 
ferroelectric distortion. In ABO3 compounds, ferroelec- 
tricty is usually related to the O 2p - B d hybridization 
[9] and Fillipetti and Hill |4] argued that because of the 
partial d-shell occupancy of Mn in CaMnOs, Mn-ion is 
unable to accept a transfer of charge from the neigh- 
boring oxygen so that there is a strong resistance of the 
filled d-shells with spherical symmetry to move off-center. 
The fact that ferroelectricity requires d^-ness while mag- 
netism requires partial d-state occupancy is often evoked 
to explain the scarcity of compounds combining both 
properties [2] |3l H] |10]. It also suggests that the search 
for multiferroic ABO3 perovskites should be oriented to- 
ward materials in which magnetism and ferroelectricity 
are driven independently by the A and B cations (3| 111). 

Khomskii et. al |12] further argued that since Mn ion 
is in the Mn^+ state, the empty e^ orbitals could never- 
theless take part in the charge transfer process. However, 



the authors suggested that the exactly half-filled or- 
bitals might restrict the bonding between empty eg or- 
bitals and oxygen p-orbitals by imposing Hund's rule on 
the eg states. More recently, a first-principles study of 
the othorhombic phase of CaMnOs |13] revealed that the 
static dielectric constant in that phase is very large and 
comparable to that of isostructural CaTiOs [S]. Since 
the latter is an incipient ferroelectric and its large dielec- 
tric response is coming from the softening of the highly 
polar ferroelectric mode, this is questioning the eventual 
tendency of CaMnOs to similarly exhibit tendency to fer- 
roelectricity. 

In this letter, we first carefully reinvestigate from first- 
principles the structural instabilities of cubic CaMnOs. 
We point out that, contrary to common believe, CaMnOs 
do develop a weak ferroelectric (FE) instability at its 
equilibrium volume. However, like in CaTiOs, the latter 
is hidden by a much deeper antiferrodistortive (AFD) 
instability which is responsible for the orthorhombic 
ground-state. FE and AFD instabilities are however in 
competition in CaMnOs and we then explore the pos- 
sibility of making it multiferroic using either strain or 
chemical engineering. 

Our first-principles calculations have been performed 
in the framework of Density Functional Theory (DFT) 
using the ABINIT package [15]. In all our calculations 
we imposed a G-type antiferromagnetic order and worked 
within the Generalized Gradient Approximation (GGA), 
using the functional recently proposed by Wu and Co- 
hen |16) . which is known to be very accurate in predict- 
ing the unit-cell volume of solids. This choice is par- 
ticularly important in the present study since ferroelec- 
tricity is strongly sensitive to the volume and accurate 
volume estimate is mandatory to make trustable predic- 
tions. We used optimized pseudopotentialsjl^ generated 
with OPIUM code|T8], treating 3s,3p,3d,4s states as va- 
lence states for Mn, 3s,3p,4s states for Ca and 2s, 2p as 
valence states for O. The wave-function was expanded in 
plane-waves, upto a kinetic energy cutoff of 55 Ha. Inte- 
grals over the Brillouin zone were replaced by sums on a 
6x6x6 mesh of special fc-points. The phonons, dielectric 
tensors and Born effective charges have been calculated 
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in the framework of density functional perturbation the- 
ory |T9] . 

First, we reinvestigated the structural, electronic and 
dielectric properties of CaMnOs in the cubic perovskite 
structure. In this high symmetry phase, the atomic po- 
sitions are fixed by symmetry so that the only structural 
degree of freedom is the lattice parameter. Our calcu- 
lations provide a relaxed lattice constant = 3.74 A, 
which is in excellent agreement with the experimental 
estimate of 3.73 A[20], illustrating the accuracy of the 
Wu and Cohen functional for magnetic compounds. 
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mized volume we surprisingly identified a small FE insta- 
bility related to a unstable polar mode with an imaginary 
frequency lofe = 13z cm~^. However, we also observed, 
within the same structure, a much larger AFD instability 
associated to a non-polar oxygen rotational mode with 
an imaginary frequency luafd — 213z cm~^. It is in fact 
the condensation of this latter mode with different ampli- 
tudes along the different directions that generates the ob- 
served orthorhombic ground-state. The case of CaMnOa 
is therefore very similar to that of SrTiOa or CaTiOs: it 
develops both FE and AFD instabilities and the absence 
of a ferroelectric ground-state is not directly related to 
its magnetic character but to a much larger tendency to 
oxygen rotations. Both types of instability are in compe- 
tition and we explicitely checked that once the dominant 
AFD distrortion is condensed, the small FE instability is 
suppressed by anharmonic effects. 



TABLE I: Lattice constant, energy gap, Born effective charges 
and optical dielectric constant of cubic CaMnOa at three dif- 
ferent volumes. 0± and Oy refer respectively to O displace- 
ments parallel and perpendicular to the Mn-0 direction. 

Within GGA, this relaxed structure is correctly pre- 
dicted to be insulating. Although significantly smaller 
than the experimental value (3.1 eV[21]) because of 
the well-known DFT bandgap problem, our calculated 
bandgap of 0.54 eV nevertheless agrees with previous 
LDA calculations jH [22] ■ On top of calculations at the 
optimized lattice parameter, we also considered struc- 
tures at a_2% — 0.98 Cq — 3.67 A and a^2% — 1-02 — 
3.81 A. The results are summarized in Table I. We ob- 
serve that the electronic bandgap decreases when increas- 
ing the volume, which is expected from the evolution of 
the crystal field in a compound where the gap arises from 
the splitting between and eg orbitals. CaMnOa nev- 
ertheless remains an insulator at the three volumes, even 
at the GGA level. 

The Born effective charges (Z*) and optical dielectric 
constants (eoo) are also reported in Table I and com- 
pare well with previous calculations |H [22] ■ We notice 
that the Born effective charges present features very sim- 
ilar to those of ferroelectric perovskites [25, with Z^,j„ 
and Zq|| (for O displacement parallel to the Mn-0 direc- 
tion) being anomalously large. In prototypical ferroelec- 
tric perovskites like BaTiOa, these anomalous Z* were 
shown to be a key feature to produce a giant destabi- 
lizing dipolar interaction at the origin of the ferroelectric 
instability [21]. Although the dielectric constant is larger 
in CaMnOa as properly noticed in Ref. [4] and will so bet- 
ter screen the Coulomb interaction, it remains surprizing 
that CaMnOs does not exhibit any tendency to ferroelec- 
tricity in spite of its giant Z*. 

We therefore computed the phonon frequencies within 
the cubic perovskite structure in order to identify the 
structural instabilities inherent to CaMnOa. At the opti- 
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FIG. 1: (Color online) Top: Evolution of the square of 
the frequency of FE mode (red circles) and AFD mode (blue 
squares) of cubic CaMnOa with the lattice constant. Below: 
Decomposition of the square of the FE frequency (red circles) 
into its short-range (blue diamonds) and long-range contribu- 
tions (blue squares) as explained in the text. 



Performing a similar characterization at different vol- 
umes, we observe in Fig. 1 that the AFD instability 
is rather independent of the lattice constant. At the op- 
posite, the FE instability is strongly sensitive to volume, 
a result totally similar to what was previously reported 
for titanates [1]. This might explain why no FE instabil- 
ity was previously reported from LDA calculations which 
typically underestimates the volume by 1-2 %. 

Following the original idea of Cochran[25], the fer- 
roelectric instability in perovskite oxides is usually ex- 
plained in terms of a close competition between short- 
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range forces, which tend to favor the high symmetry 
phase, and long-range Coulomb interactions which tend 
to destabilize the structure. In order to clarify wether 
this idea still apply to CaMnOs or wether its ferroelectric 
instability has a different origin, we performed a decom- 
position of the FE frequency in terms short-range (SR) 
and long-range (LR) Coulomb contribution following a 
scheme previously applied to BaTiOa in Ref. p4] : 

{lFE\Dtot\lFE) = {ife\Dlr\ife) + {ife\Dsr\ife) 

V ' " V ' " V ' 

'^'fe ^lr ^Ir 

where Dtot refers to the full dynamical matrix at F, Dlb. 
and DsR to its LR and SR contributions and |7f_e) is 
the dynamical matrix eigenvector of the FE mode. The 
results are summarized in Figure [2] and are comparable 
to those previously reported for BaTiOs : the small FE 
instability is a consequence of a delicate balance between 
SR and LR forces and arises from a giant negative LR 
contribution directly related to the anomalously large Z* 
values. We explicitely checked that, considering them at 
their own optimized volume, CaMnOa develops a DD 
contribution {lodd = 378. 4i cm~^) even larger to that 
of CaTiOa {lodd = 345i cm~^) and that in spite of its 
larger dielectric constant. The smaller FE instability of 
CaMnOs is related to stronger SR forces associated to 
its smaller equilibrium volume {u)sR — 378.2 cm~^ in 
CaMnOs versus losr — 313 cm~^ in CaTiOa) . 

The natural tendency of CaMnOs to develop a fer- 
roelectric distortion being established, it might now be 
questioned which cation is driving the FE instability. To 
clarify this issue, the involvement of each type of atom 
within the ferroelectric soft-mode eigendisplacement vec- 
tor is reported in Table II [26]. Inspection of this Ta- 
ble points out that the progressive softening of the FE 
mode when increasing the volume is closely related to 
the concomitant increase of the Mn atom motion. This 
involvement of the Mn atom which has the largest Z*, 
combined to a lower extent with the small increase of Z* 
itself, produce a significant increase of the mode effective 
charge and of the related destabilizing DD interaction. 
This clearly establishes the Mn driven character of the 
ferroelectric instability. CaMnOs provides therefore a 
first prototypical example of magnetic perovskite ABO3 
compound in which ferroelectricity and magnetism are 
driven by the same cation, illustrating the fact that mag- 
netism and ferroelectricity are not necessarily exclusive. 

The absence of FE phase transition in CaMnOs being 
related to the competition with the AFD instability, we 
now investigate to which extent it would be possible to 
favor ferroelectricity and so, eventually, make CaMnOa 
multiferroic using either strain or chemical engineering. 

From Fig.l, it appears that FE and AFD instabilities 
exhibit a very different sensitivity to strain, a fact that 
can be potentially exploited to favor the ferroelectric dis- 
tortion. Since applying a negative isotropic pressure is 



ABO3 acell l^AFD <^FE Zp^ 2??^ 2?7b Qrjo 

A cm~^ cm~^ |e| 

CaMnOs 3.67 223i 146 7.51 0.59 0.00 0.41 

3.74 213i 13i 10.22 0.34 0.05 0.61 

3.81 205i 189i 13.48 0.08 0.26 0.66 

SrMnOs 3.81 55i 148i 14.65 0.00 0.19 0.81 

Cao.sBao.sMnOa 3.83 150i 216i 15.28 0.00 0.21 0.79 

BaMnOa 3.92 140 335i 19.24 0.00 0.22 0.78 



TABLE II: Lattice constant (flceii), frequencies of the AFD 
{ujafd) and FE [ujfe] phonon modes, mode effective charge 
of the FE mode (Zpp) and involvement of each type of 
atom to the FE mode eigendisplacement vector |26] for dif- 
ferent compounds in the cubic perovskite structure. For 
Cao.sBao.sMnOa, results correspond to an ordered supercell 
in witch Ca and Ba alternate along the [111] direction. 



not easily achievable in practice, we rather explore the 
possibility to make CaMnOa ferroelectric under epitax- 
ial strain as it could be achieved in epitaxial thin films. 
In ferroelectric oxides, it is well-known |27) that large 
enough compressive epitaxial strains typically favor a 
ferroelectric c-phase (with the polarization out-of-plane) 
while tensile strains favor a ferroelectric a-phase (with 
polarization in-plane) but no systematic data are avail- 
able for the impact of epitaxial strain on the AFD distor- 
tion. Exploring first the effect of compressive epitaxial 
strains, we observed that the ground-state of CaMnOa 
remains AFD distorted and unfortunately paraelectric 
for strains as large as —4%, roughly the largest epitax- 
ial strain experimentally achievable. At the opposite, for 
tensile strains, it appeared however that a polar ground- 
state is obtained for strains larger than about -1-2%. Be- 
yond this critical value, the system preserves its oxygen 
rotations by aquires an additional in-plane polarization 
estimated from a Berry phase calculation to 4 /iC.cm"^ 
at an in-plane lattice constant of 3.82 A. This offers the 
theoretical demonstration that CaMnOa can be made fer- 
roelectric by strain engineering. 

Another possible way to increase the volume and so 
provide more space for the FE distortion of the Mn atom 
is to perform atomic substitution and to replace Ca by 
a larger ion. To support this idea we compare in Figure 
2 and Table II the behavior of CaMnOa, SrMnOa and 
BaMnOa, each of them in the same prototypical cubic 
perovskite structure at their own equilibrium volume. It 
clearly appears that increasing the size of the A cation 
strongly favor the FE distortion that is the only remain- 
ing instability in BaMnOa. Going further we see that 
the atomic substitution goes beyond a simple volume ef- 
fect since together with favoring ferroelectricity, bigger 
A cations also suppress the possibility of oxygen rota- 
tion. From that viewpoint, atomic substitution of Ca by 
a bigger atom appears even more promizing than epitax- 
ial strain to produce a ferroelectric ground-state. 

To support this idea, we performed calculations on al- 
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loy Cao.sBao.sMnOs. Working first within the virtual 
crystal approximation (VGA) we obtained a relaxed 
lattice constant ao = 3.84 A and observed a strong FE 
instability {ujpE — 224z cm~^) while the AFD instability 
was totally suppressed {luafd = 110 cm~^). Consider- 
ing then an ordered supercell in which Ca and Ba ions 
alternates along the [111] direction, the relaxed lattice 
constant (ao = 3.83 A) and the ferroelectric instability 
{copE = 215i cm"^) were very similar to VGA but the 
AFD mode was now unstable {ojafd = 150z cm~^). This 
points out that the amplitude of the AFD mode might be 
very sensitive to the atomic order but both approaches 
predict a dominant ferroelectric instability. This con- 
firms that, eventually combined with strain engineering, 
partial substitution of Ca by Ba will favor ferroelectric- 
ity and might be also a promizing way to make CaMnOa 
multiferroic. 




3.70 3.75 3.80 3.85 3.90 3.95 

Lattice constant (A) 

FIG. 2: (Color online) Square of the frequency of FE (red 
circles) and AFD (blue squares) modes of CaMnOs (CMO, 
ao =3.74 A), SrMnOa (SMO, ao =3.81 A) and BaMnOa 
(BMO, ao =3.92 A) at their own equilibrium volume. 



In conclusion, we have investigated from first-principles 
the lattice dynamics of cubic CaMnOa. We pointed out 
that, together with a strong AFD instability responsible 
for its orthorhombic ground-state, CaMnOa also surpris- 
ingly exhibits a weak FE instability at its equilibrium vol- 
ume. Going further, we highlighted that both strain and 
chemical engineerings can be used to favor the ferroelec- 
tric distortion and make CaMnOa multiferroic. Finally, 
we also showed that the FE instability of CaMnOa is 
Mn-dominated, illustrating that ferroelectricity and mag- 
netism are not necessarily exclusive and can be driven by 
the same cation. All this suggests alternative strategies 
in the search of multiferroic ABOa compounds and we 
hope that our findings will motivate further experimen- 
tal investigations. 
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